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Abstract 

Foliate systems are those which preserve some (possibly singular) foliation of phase 
space, such as systems with integrals, systems with continuous symmetries, and skew 
product systems. We study numerical integrators which also preserve the foliation. 
The case in which the foliation is given by the orbits of an action of a Lie group 
has a particularly nice structure, which we study in detail, giving conditions under 
which all foliate vector fields can be written as the sum of a vector field tangent 
to the orbits and a vector field invariant under the group action. This allows the 
application of many techniques of geometric integration, including splitting methods 
and Lie group integrators. 

1 Introduction 



In the early works of Feng Kang [8] , geometric integration was taken to be the 
approximation of flows by elements of certain subgroups of Diff{M) (the group 
of diffeomorphisms of the phase space M) — the groups of symplectic, volume- 
preserving, or contact diffeomorphisms, for example. This point of view was 
developed further in [14] using the Cartan classification of diffeomorphism 
groups and is continued here by considering the so-called nonprimitive groups, 
those that leave a foliation of M invariant. (A preliminary announcement of 
some of our results appears in [15].) That is, we ask: given a vector field whose 
flow preserves a given foliation of M, how can we construct integrators with the 
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same property? Many other interesting (and difficult!) questions can be asked 
aboTit this class of systems. How can the existence of an invariant foliation be 
detected? What are the consequences for the dynamics of the system? Even 
regarding the construction of integrators, the class of all foliations seems to 
be too large to admit a useful theory, and we are led (following the example 
of Lie group integrators [9,21]) to consider foliations defined by the action of 
a Lie group. We introduce these with an example. 

Example 1 Let M = and consider the vector field 

X — xy + x{l — x^ — y'^), y — —x^ + y{l — x^ — y^) (1) 

In polar coordinates, this becomes 

r = r(l — r^), 9 = —rcos9, 

showing that the foliation into circles r — const, is invariant under the flow. 
(In fact, this foliation is singular, because the leaf through the origin, a single 
point, has less than maximal dimension.) A one-step integrator preserves this 
foliation (i.e., is 'foliate' if the final value of r is independent of the initial 
value of 9. Of course this is easy to obtain in polar coordinates, but we shall 
see that no standard integrator in cartesian coordinates is foliate. The leaves 
of this foliation are the group orbits of the standard action of 5*0(2) on 
(see Figure 1, top). 

We show in section 3 that in many cases, vector fields whose flows preserve such 
a foliation have a particularly simple structure, namely, they can be written 
as a sum of two vector fields, one tangent to the leaves of the foliation (being 
the orbits of a group action), and one invariant under the group action. This 
representation allows the construction of various types of foliation-preserving 
integrators. 

Definition 1 [20,22] Let M he a manifold of dimension m. A singular folia- 
tion F of M is a partition of M into connected immersed submanifolds ( the 
"leaves"), such that the vector fields on M tangent to the leaves are transi- 
tive on each leaf [22]. F is regular if each leaf has the same dimension. F 
has codimension q if the maximum dimension of the leaves of F is m — q. 
A dijfeomorphism of M is foliate with respect to F if it leaves the foliation 
invariant, i.e., if it maps leaves to leaves. A vector field on M is foliate if its 
flow is foliate. The space of smooth vector fields foliate with respect to F is 
denoted Xf- A (one-step) integrator is foliate with respect to F if the dif- 
feomorphisms of M corresponding to each time step are foliate. The space of 
smooth vector fields tangent to the leaves of F is denoted Xtan- The space of 
leaves ( denoted M/ F) is obtained by identifying the points in each leaf together 
with the quotient topology. 
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Fig. 1. Three foliate vector fields. Top: a general foliate vector field, r = r(l — r^), 
9 = —rcosO (Eq. (1)). Middle: a system with an integral, f = 0, 9 = —rcos9. 
Bottom: a system with a continuous symmetry, r = r(l — r^), 6 = —(1 + r^/5). All 
three flows map circles to circles. The dots mark times 0, 0.5, and 1. 



3 



Theorem 1 [20] Xp and Xtan form Lie algebras. Xtan is an ideal in Xp- A 
vector field X is foliate with respect to F if and only if [X, Y] e Xtan for all 

Y G ^tan- 

Theorem 2 [20] Let M and N be manifolds of dimension m and n, respec- 
tively. Let I : M ^ N be a smooth surjection. (If I is not onto, we replace N 

by I{M).) Then I defines a foliation F whose leaves are given by the connected 
components of I^^{y) for each y G I{M). If I is a submersion, i.e. if dl has 
constant rank n, then F is a regular foliation of codimension n. In this case 
the space of leaves M/F is diffeomorphic to N. 

Such a foliation is called simple. Given a vector field, one can search for simple 
foliations it preserves by looking for functions / such that / = /(/). These 
functions / might be conserved quantities in special cases. For example, the 
conserved momentum function J of a Hamiltonian system can evolve under 
J — cJ under the addition of linear friction. Then the perturbed system 
preserves the simple foliation defined by J [18]. 

For a regular simple foliation, a foliate vector field X drops to a vector field 
on the space of leaves M/F, which (following the tradition in geometric me- 
chanics, [13]) we call the reduced system. If ^ : [0,T] — > M/F is an orbit of 
the reduced system, the problem of finding an orbit of X is called the recon- 
struction problem. If all leaves are diffeomorphic, then reconstruction involves 
integrating a nonautonomous vector field on a fixed leaf. (In Example 1, the 
reduced system is r = r(l — r^); for any solution r{t) of this equation, the 
reconstruction system is ^ = ~^(^) cos^, a nonautonomous ODE on the leaf 
S\) 

In the application to integrators, we do not usually want to construct the 
reduced system explicitly since the original phase space M is usually linear 
and easier to work in. We would like integrators on M which preserve the 
foliation automatically. 

Example 2 A system with k first integrals / : M — >■ is foliate with respect 
to the level sets of the functions /. Each leaf is in fact fixed by the fiow. For 
this reason we choose the symbol / in Theorem 2 to suggest that simple foliate 
systems generalize systems with first integrals. (See Figure 1, middle.) 

Example 3 A system with a symmetry is foliate with respect to the orbits of 

the symmetry. That is, let X admit the Lie group action A:GxM— >Masa 
symmetry, so that its fiow 0^ is G-cquivariant. Then \{g,(j)t{x)) = ([>t{\{g , x)) , 
i.e., the foliation with leaves given by the group orbits {X{g,x) : g E G} is 
invariant. In this case the reconstruction problem on G is easier to solve than 
in the general case, because it is G-invariant. (See Figure 1, bottom.) 
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Example 4 The Lorenz system is given by 



X — ay — ax, 

y — —y — xz — rx, 

z — xy — bx. 

If 6 = 2(7, the system is fohate with leaves x"^ — 2az = const., for 

^(x^ - 2a z) = -2a{x'^ - 2a z). 

(A fohate integrator can be constructed as foUows. We spht into 

Xi. X = ay, y — —xz — rx, z = xy, 
X2. X — —ax, y — —y, z — —2ax. 

Xi is tangent to the fohation and may be integrated using the midpoint rule, 
which preserves the quadratic function — 2az; X2 is foliate but linear, and 
can be solved exactly. Composition then yields a foliate integrator.) 

Example 5 A special case of the foliations defined by submersions is given 
by M — N X L and / is projection onto N. Each leaf is then diffeomorphic 
to L. (This is locally true of any simple foliation in a neighborhood in which 
the leaves have constant dimension.) In coordinates x on N and y on L, any 
foliate vector field can be written in coordinates as 

X = f{x) 



and any tangent vector field as 



x = 

y = g{x,y). 



These foliate vector fields arc also known as skew product systems, introduced 
by Anzai [1], and studied today in ergodic theory [3] and complex dynamics 
[10]. A special case of skew product systems is given by iV = M" and L = M™-"", 
the standard foliation on W". 



Example 6 The extension of a nonautonomous vector field on M to an au- 
tonomous vector field on M x M preserves the foliation defined by t = const. 
Most integrators are foliate and indeed, solve the reduced system t — 1 exactly. 

In a foliate system, one can obtain some information about part of the system 
(namely, the current leaf) for all time without even knowing the full initial 
condition. Surely this puts strong dynamical constraints on the whole system. 
Nevertheless, the only properties we can point to depend on the reduced or 



5 



reconstruction systems being simple in some way: having small dimension or 
simple (e.g. linear) dynamics. 

Example 7 Burns and Wilkinson [3] study skew product systems of the form 
{x,g) I— s> {f{x),(f)(x)g) on M X G. where G is a compact Lie group, : M ^ 
G, and / is measure preserving. Here the foliation is simple as in Example 
5 and the systems have extra structure: the reduced system x i-^ f{x) is 
measure preserving, the reconstruction system g i— > 4>{x)g is G-equivariant, 
and the whole system is measure preserving (with respect to the product of 
the measure on M with Haar measure on G). This class is a special case of 
the class of measure preserving systems with a symmetry. The extra structure 
allows them to prove that the ergo die systems are open and dense in this class. 

Example 8 Consider the system on with a codimension 1 foliation, 

X = f{x), y = g{x, y,z), z = h{x, y, z). 

The only possible c<j-limit set of the reduced system i = /(x) is a point, which 
suggests that in this case the cj-limit set of the whole system is either a point or 
a circle (periodic orbit). (If f{x) > 0, however, so that the reduced system has 
no a;-limit set, chaos is possible in the full system.) Similarly, for the system 
on with a codimension 2 foliation, 

X = f{x, y), y = g{x, y), z = h{x, y, z). 

the possible c<j-limit sets of the reduced system in (x, y) are a points and circles, 
which suggests that in these cases the u;-limit set of the full system is a point, a 
circle, or a 2-torus. In both examples the existence of the foliation influences 
the possible dynamics. 



2 Integrators for simple foliations 

In this section we consider whether standard integrators can be foliate, and 
adapt geometric integrators for systems with integrals to the foliate case. 

The next Theorem generalizes the fact that Runge-Kutta methods preserve 
arbitrary linear integrals. 

Theorem 3 Let M = and let F be a linear foliation, i.e., a simple fo- 
liation defined by the linear function I : R™ — > R'^. Then any Runge-Kutta 
method is foliate. 

PROOF. Runge-Kutta methods are linearly covariant [17], hence we can 
apply a linear change of variables to bring I into the form I{xi, . . . ,Xm) — 
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(xi, . . . ,Xk). Runge-Kutta methods are closed under restriction to the closed 
subsystem with coordinates Xi, . . . ,Xk [2], hence the final values oi Xi, . . . ,Xk 
depend only on their initial values, i.e., the method is foliate. 



We now give two methods for constructing foliate integrators, each based on 
a popular integral-preserving method. 

Projection methods. Let the leaves of the foliation be the level sets of the 
functions 7 : M — > R'^ and write the reduced system as 7 = h{I). Let x„ be the 
initial condition. Then the following algorithm provides a foliate integrator. 

Step 1. Calculate /„ = /(a;„) and integrate the reduced system for one time 

step, giving which depends only on /„. 
Step 2. Integrate the full vector field for one time step, giving x. 
Step 3. Calculate Xn+i by projecting x onto the desired leaf I~^{In+i), e.g. by 

orthogonal projection. 

Discrete gradient methods. [16] For simplicity we present the method for a 
single function I, i.e., for a codimension 1 foliation, and use a Euclidean metric. 

It is necessary to first split the vector field into its components tangent and 
orthogonal to the leaf. Since I = X.VJ = h^I), the component of X in the 
direction V/ is (/i(/)/|V/p)V/, while the component of X tangent to the leaf 
can be written in the form y4(a;)V/ for some antisymmetric matrix A{x) [16]. 
That is, the full system can be written 

Extend A and h to functions A and h of pairs of points satisfying A antisym- 
metric, A{x, x) = A{x), h{I, I) = h{I), and let V be a discrete gradient. Then 
the discrete-gradient discretization 

obeys 

1{X ) - I{x) = Th{I{x),I{x')), 

SO the integrator is foliate. Here r is the time step. 

One popular integral-preserving method does not generalize to foliations. The 

symplectic Runge-Kutta methods preserve arbitrary quadratic first integrals 
(and arbitrary linear symmetries) but do not preserve arbitrary quadratic fo- 
liations. To see this, consider the midpoint rule applied to a system preserving 
the foliation r = const, as in Example 1. Let the method he x ^ x' and write 
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X — {x + x')/2. The method is x' — x + tX{x). We have 

r'^ - = {x' + xYix' -x) = 2tx^X{x). 
Since the vector field X is fohate, 

4--r^ = = x^X(x) 
dt2 ^ ^ 

is a function of r only. Therefore, r'^ — is a function of \x\'^ only. However, 
|xp is not in general a function of r. For example, consider the system 

x^ -y^ + x, y^xy + y (2) 

for which r = r, ^ = rsin^. Applying the midpoint rule and expanding in a 
Taylor series gives 

r'' = r'(l + 2r + 2t^ + t^{3 - y^)/2) + 0(r^) 

which is not a function of r only. 

To obtain foliate integrators which are 'intrinsic' in the sense that they do not 
involve constructing the reduced system explicitly, we have to consider special 
cases. For example, linear foliations can be preserved in a sense automatically, 
and for systems with linear symmetries, the foliation defined by the orbits of 
the symmetry is preserved by symmetry-preserving integrators. We therefore 
consider in the next section foliations defined by the orbits of a Lie group 
action. Not only is the foliation simple, but the foliate vector fields themselves 
have a nice structure which allows the construction of foliate integrators. 



3 Lie group foliations 

Let G be a Lie group and A: G x M — > M be an action of G on M. (We write the 
action as \{g, x) — XgX — g • x — gx, as needed.) This group action generates 
a (possibly singular) foliation whose leaves arc the group orbits \{G,x). A 
vector field preserving this foliation is said to be G-foliate. Let q be the Lie 
algebra of G and let Qm be the distribution tangent to the leaves, with each 
i E g associated to a vector field on M, i.e., 

Qm{x) = T,{X{G,x)) = {^ix):^ e g}. 

(Recall that p-dimensional distribution on M is given by associating, to each 
point X in M, a p-dimensional subspace of the tangent space Tj.M [20].) 
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As with any foliation, any vector field tangent to the leaves is foliate. However, 
in the case of a Lie group foliation we have another natural class of foliate 
vector fields, namely those invariant under the action. For, from Theorem 1, 
X is foliate iff [X, Qm] C 0m; but if X is G-invariant, then [X, Qm] = 0. Thus, 
the vector field Xtan + Xinv is G-fohate, where Xtan is tangent to the leaves and 
Xinv is G-invariant. To characterize the G-foliatc vector fields we now want to 
turn this around and ask: When can all foliate vector fields be decomposed 
in this way? We shall see that this decomposition makes it easy to construct 
foliate vector fields exphcitly and also to construct foliate integrators. 

Our results can be summarized as follows. Let X be a foliate vector field. 
First, if G acts by isometrics, then X = Xtan + X;j„ where Xjnv is not only G- 
invariant but also perpendicular to the leaves (Theorem 5). Second, if has 
a G-invariant complement if, then X — X^^^ + Xjnv with Xjnv G H (Theorem 
6). This arises, for example, if the action is free and proper (Theorem 7). The 
decomposition is also true locally in a neighborhood of an orbit which admits 
a certain slice (Theorem 8). Finally, we give a counterexample to show that 
not all actions admit such a decomposition (Example 9). 

Recall that a proper action is one for which GxM ^ MxM, {g,x) i— > {g.x, x) 
is proper, i.e., the inverse image of any compact set is compact. An action of a 
compact group must be proper. A proper action must have compact isotropy 
groups at all points of M. 

Theorem 4/7/ Let G he a Lie group acting properly on the manifold M. 
Then M has a G-invariant Riemannian metric. 

Theorem 5 Let (M, (, )) he a Riemannian manifold. Suppose we have a smooth 
group action X:GxM^M acting by isometrics on the metric (, ). Consider 
the foliation hy the group action where the leaf through x & M is simply the 
orhit X{G,x). Then the vector field X is G-foliate if and only if the unique 
metric decomposition of X into tangential and perpendicular components 



satisfies, for all g G G, (t)*gX = X . In other words, X is G-foliate if and 
only if its component perpendicular to the leaves is G-invariant. 



PROOF. Recall that Xtan is the space of smooth vector fields tangent to the 
foliation. In this case we can write 



X = xii +X 



■± 



(3) 



(4) 
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Prom Theorem 1 we have that X is fohate if and only if for all Y e Xtan, 

[x,F] = [xii+x^,r]eXtan. (5) 

However, G Xtan by definition, so X is foliate if and only if for all Y e Xtan, 

e^tan. (6) 

Now, suppose X is foliate. The condition above requires, in particular, that for 
all ^ G 0, [^Mi G Xtan- Now, by definition of X^ we have that (x^, ^m) = 
0. Fix X E M and t] E Q. Let (j)} denote the flow of r]M which is given by 
4>t{x) = \{exptri,x). We then have 

x^{rt(^),^Mm^)))^o. (7) 

Since (pf leaves the metric invariant for all t, we have 

[Tcl>l,X\^'l{x)),T^lM'Pn^))) = (8) 



so that 

^' (T0!L,X^(0?(x)),T0l,eM(0?(a;)) 
= {L^^X^{x),^m{x)) + (X^{x), L^^Uix)). 



dt 



(9) 



Since the action is on the left we have L^^^^^m = [vm,Cm] — —[VtC]m, so that 
the second term in Equation (9) is 0. It follows that 

L,^X^,eM) = (10) 



which shows that Lj^^^X-^ is perpendicular to the leaves. However, by the 
assumption that X is foliate we demand that [X-^,r]M] = L^^X-^ G Xtan 
so that we must have L^^^X-^ — 0, proving G-invariance of X-^. Conversely, 
suppose X"*- is G- invariant. It suffices to show that X-^ is foliate with respect 
to the G orbits. But this is immediate since the fiow of X-^ is G-equivariant 
and therefore, if 0t denotes the flow of X"*-, 4>t^{g, x) = \{g, 4>t{x)) which shows 
that (pt maps the leaf through x to the leaf through (f)t{x). 

Theorem 6 Let \: G ^ M be a smooth action of the Lie group G on the 
manifold M. Suppose the tangent bundle TM admits a G-invariant splitting 
TM = qm + H. Let the corresponding decomposition of a vector field X be 
X — Y + Z . Then X is G-foliate iff Z is G-invariant. 
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PROOF. The 'if part is immediate. Conversely, let X be a G-foliate vector 
field. We have to show that Z is G-invariant, i.e., that [^Mi -Z^] = for all ^ G 0. 
Since X is fohate, [Cm,^] = [Cm, 5^ + ^] G 0m; but [Ca^,F] G 0m since Y is 
tangent to the leaves by definition. Therefore [^m, ^] G 0m- Fix x E M and 
g & G. Since if is G-invariant, we have 



Letting = exp{t^M) and differentiating with respect to t at t — gives 



Since the two subspaces Qm{x) and ifj. are complementary, we must have 
[^M, Z] = 0, that is, Z is G-invariant. 

The next theorem follows either from Theorem 5 or from Theorem 6. How- 
ever, we also give a direct proof which constructs the decomposition explicitly 
relative to the choice of a connection. Recall that a free action is one in which 
all isotropy groups are trivial, Gx — {e} for all x E M. 

Theorem 7 Let X: G x M ^ M be a smooth, free, and proper action of the 
Lie group G on the manifold M. Let X be a G-foliate vector field. Then there 
exists a vector field Xtan G 3£tan and a G-invariant vector field Xi^^ such that 

X — ^tan + ^inv 

PROOF. Because the action is free and proper, M — > M/G is a principal 
G-bundlc. There is a 1-1 correspondence between the space of connections 
and G-invariant complements to 0m- Fix a connection A and let H^. = kcr^^-. 
Furthermore, there is a unique horizontal lift map hor: T(M/G) — > TM. It 
then follows that horvr^X is a G-invariant vector field taking values in H. 
Thus Xinv = hoTTT^X and Xtan = — -'^inv provides the desired decomposition 
of X. Xtan is tangent to the foliation because ^(-'^tan) = A{X — hor7r*X) = 
A{X) - A{horn,X) = A{X). 

Definition 2 [7] A slice at Xq E M of a smooth Lie group action X:GxM^ 
M is a submanifold S of M through Xq such that 

(i) T^gM = Qm(xo) © T.^oS (^'^d for all x e S, T^M = Qm{x) + T^S; 
(a) S is GxQ-invariant; and 

(Hi) if X & S , g & G, and \{g,x) G S, then g G G^^- 



Tg.xXg-i{Z{g ■ x)) G Hx. 



[Cm,Z]{x) = L^^^Z{x) 




'exp(t^)'^(^) 



Te^pit^)-xKxpi-t^yx{Z{ex-p{t^) -x)) G H^. 
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Note that G ■ S is an open, G-invariant neighborhood of the orbit G ■ xq, 
and every orbit in this neighborhood intersects in a unique G^(,-orbit of S. 
Every proper action admits a shce at every point. The dimension of an orbit 
which intersects S must be greater than or equal to the dimension of the orbit 
through X. 

Unfortunately, it appears that existence of a slice is not sufficient to guarantee 
the existence of a G-invariant splitting of T{G ■ S). The following theorem 
requires two extra assumptions. However, the first is a requirement on the 
group itself (not on the action), while the second only concerns the action 
of Gxo on S, reducing the dimensionality. One important special case (to be 
illustrated in Example 10) is when Gx = Gx^ for all x in S, in which the second 
assumption is automatically satisfied with H = TS. The following theorem 
gives sufficient conditions which, together with Theorem 6, guarantee that on 
G-S every G-foliate vector field can be decomposed into tangent and invariant 
components. 

Theorem 8 Let X:Gx M ^ M be a smooth action of the Lie group G on the 
manifold M which admits a slice S through the point xq & M. Furthermore, 
assume 

(i) Q admits a G^Q-invariant splitting Q = Qxo®^ where the group GxQ acts on 
by adjoint action. 

(a) There exists a G^Q-invariant splitting TS — {Qxo)m ® H of the tangent 
bundle of S. 

Let U be the open G-invariant neighborhood G • S of the orbit through xq. Then 
the tangent bundle TU admits a G-invariant splitting 



T^M = g-H,® {g-i)M («) © {Q9-xo)m («) 



fill 



where u G U is given by u = gs, the first distribution on the right hand side is 
transverse to the group orbits and the other two form a decomposition of the 
group directions into directions transverse and tangent to the slice respectively. 



PROOF. We first check that each distribution in Eq. (11) is G-invariant. 
Since the action satisfies Gg.xo = gGx^g'^, we have Qg.x^ = ^dgQxo ='■ 9 ■ Qxo- 
Therefore 

Next, since g ■ {C)m ~ {9 ' Om holds for any group action, we have 

TXg:{t)^{s)c^{g-t)^{gs). 
It follows that {g ■ t)^ {gs) fl {g ■ Qxo)m (9-^) ^ ^- Next, the first distribution is 
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G-invariant by construction, and since TXg-. Ts{G ■ s) ~ Tgs{G ■ gs), we must 
have^-i/,nr^,(G'-^s) = 0. 

Finally, we must check that the splitting is well defined. Suppose u — gs — 
giSi. Then without loss of generality we can take si = s and gi = gh with 
h G GxQ, since by the shce property, Gg C G^q- (Each orbit in U intersects 5* 
in a unique G^^g-orbit.) The sphtting is then weU-defined if TXh maps each of 
the three distributions isomorphically into themselves. We clearly have that 
TXh-.T.S ~ ThsS and that TXh:Ts{G^, ■ s) ~ Ths{G^, ■ s), since h E G^, . So 
what we need is that TXh maps Hg into Hhs, which is true by assumption (ii), 
and that /i • t = t for all /i e G^o , which is true by assumption (i) . 



We do not have necessary and sufficient conditions for the action to be such 

that the tangent and invariant vector fields span all foliate vector fields. How 
ever, the following example shows that this is not true for all group actions. 

Example 9 Let M — R^ and A: ((a, b), {x, y)) ^ {x,y -\-a-\-hx). The orbits of 

A are the lines parallel to the y-axis, so the foliate vector fields all have the form 
f{x)-^+g{x, y)^- The tangent vector fields have the form g(x, y)-^- However, 
the vector field X = c{x, y)-^+d(x, y)^ is G-invariant iff [^m, X] = [rjM, X] = 

0, where = and r]M - ^ 
gives 



x-^ are the two generators of the action. This 



dy ' dx dy 



— rl — 



0, 



implying Cy — dy — 0, and 



[Vm^X] 



d d . d 
oy ox ay 



-c^ = 0, 
ox 



implying c = 0. That is, X 



dy' 



and the tangent and invariant vector 



fields do not span the foliate vector fields. 



The flow of a foliate vector field necessarily maps leaves to leaves diffeomor- 
phically. In particular, each orbit is restricted to leaves of constant dimension. 
However, for vector fields of the form Xtan + ^inv, even more is true. 

Theorem 9 The flow of the G-foliate vector field X — Xta,n + -'^inv preserves 
the isotropy subgroup of the initial condition up to conjugacy. Specifically, 
Gx{t) = 9if')Gx{o)gii)~^ , where g{t) satisfies Eq. (13) below. 



PROOF. We write x{t) = X{g{t),h{t)) as in Eq. (12). The fiow 0t of Xin^ 
does not change the isotropy subgroup. G-equivariance means (f)t{X{g, h{0)) — 
X{g,(t)tih{0)) for all g e G, so g e Gh(o) iS g e Gh{t)- Then we have G^^t) = 
Gg{t)h{t) = g{t)Gh{t)g{t)-^ = g{t)Gh{o)g{t)-^ = g{t)Gx{o)g{t)-^, as required. 
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Note that the result is not true for all G-foliate flows. In Example 9, the 
isotropy groups are G(^x.y) = {{a,h):a = —bx}, which are all invariant under 
conjugacy, because G is abelian. Therefore G(^xi,y) is not conjugate to G{x2,y) 
for xi ^ X2- The flow of the fohate vector field x — 1, y — Q therefore maps 
leaves to leaves of non-conjugate isotropy. 



4 Integrators for Lie group foliate vector fields 

Given a Lie group foliate vector field X — X^an + -'^inv; this splitting into 
tangent and invariant pieces may not be unique, because it depends on the 
choice of G-invariant metric (or splitting of TM). In fact we have the whole 
family of splittings X = (Xtan + ^) + (-^inv — ^) for any G-invariant vector field 
Y e Xtan- Nevertheless, once such a sphtting has been found, it can be used to 
construct foliate integrators. First, one can integrate each piece separately and 
compose the results: Xtan by any integrator for vector fields on homogeneous 
spaces [21], and Xinv by any symmetry-preserving integrator, although these 
are easy to construct only when the action is linear. 

It is also possible to construct integrators in one piece. Write the solution 

x{t)^\{g{t)Mt)), m(0) = x(0), ^(0) = 1, (12) 

and differentiate to get 

X — DiX{g, m)g -\- D2X{g, m)fn. 

Here Di\{g,m):TgG T^M is the derivative of A in its first slot, and 
D2X{g,m):TmM — > T^M is the derivative of A in its second slot. Note that 
the first term on the right is tangent to the foliation. Therefore, we choose g 
and m so that the first term is Xtan and the second term is Xinv: 

g = iD^X{g, m))-' X,,^{X{g , m)) (13) 
m = (D^Xig, m))-^Xi„v(A(^, m)) = X,^,{m). (14) 

Here we arc using that the range of DiX{g,m) is Xta.n{x). If the action is not 
free then this map is not injective — many Lie algebra elements generate the 
same tangent vector Xtan(a^) — and some choice of the inverse {DiX{g,m))~^ 
must be made to specify a vector field on G. However, in the examples below 
a natural choice can be made. Thus we have extended the vector field on M 
to a vector field on G x M with extra foliate structure of its own. Integrating 
the reduced vector field Xi^y{m) by a symmetry-preserving method leaves a 
reconstruction problem on G. However, the vector field on G is now known 
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only at the time steps of the sohition of Xi,iv(?7i) which could make it difficult 
to achieve high orders. It is better to integrate the full system together as 
follows. 

Theorem 10 When M = R" and the group action is linear, the following 
algorithm provides a foliate integrator of order p for X = Xtan + -'^inv (i) Let 
x{t) = X{g{t),m{t)) , g{0) = 1, m(0) = a;(0) and choose a differential equation 
on G according to Eq. (13). (ii) Apply a Runge-Kutta method of order p to 
m — Xinv(?«), and the associated Runge-Kutta-Munthe-Kaas (RKMK) [21] 
method to the chosen differential equation on G, with stage values for m given 
by the corresponding stage values of the RK method applied to m — Xi-ay{m). 

PROOF. Since the group action is linear, Xi^^ has a linear symmetry and this 
is preserved by the Runge-Kutta method. That is, the integrator on m is G- 
equivariant and hence preserves the foliation into G-orbits. Applying RKMK 
to the differential equation on G is equivalent to applying the same RKMK 
method to the full system on G x M with respect to the group action of G x M 
on itself given by 

A((5'i, mi), (5f2, m2)) = (5^15^2, rui + 

and hence it has order p. Reconstructing the solution by acting on m{t) with 
g{t) only moves the point around on the new leaf, so the total integrator is 
foliate. 

Example 10 The matrix differential equation L = [A[L), L] motivated much 
of the work in Lie group integrators [9,4]. Let G C GL{n) be a matrix Lie 
group, let M — Q, the Lie algebra of G, and let Gi be a subgroup of G which 
acts on M by adjoint action, i.e. \{U, L) = ULU~^ for U E Gi, L E Q. Recall 
that the "isospectral manifolds" of Ql{n) are the sets of matrices similar by 
an element of GL{n). In the present example, the leaves of the foliation are 
the sets of matrices in q which are similar by an element of Gi, and hence are 
submanifolds of the isospectral manifolds. 

All vector fields tangent to the leaves can be written in the form 

where A: g — > 0i. From Theorem 5, when q admits an adjoint invariant metric 

(e.g. when it is compact), all foliate vector fields can be written in the form 

L^[A,L\+f{L), f{ULU-')^Uf{L)U-'yUeG, (15) 
for some (Gi)-adjoint- invariant function /:0 — > 0. 

It is not so easy to explicitly construct all invariant vector fields. The classical 
method requires determining a complete and independent set of differential 
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invariants of the action. Another approach which can work if the action is 
understood weU enough is to choose a reference point on each orbit, choose 
an isotropy-invariant tangent vector at each reference point, and push them 
around the leaves using the group action. 

For example, consider the case G = Gi = GL{n). The orbits are the conjugacy 
classes of matrices so a reference point on each orbit is given by the Jordan ma- 
trices. Restricting for simplicity to the diagonalizable matrices, whose isotropy 
groups are the group of nonsingular diagonal matrices, we first choose /(A) 
diagonal (so that it is isotropy invariant), and then define 

/(L) = Ug{A)U-' 

where L = UAU~^ is the diagonalization of L. This constructs all invariant 
vector fields on the diagonalizable matrices. 

Furthermore, this case provides a nice example of Theorem 8 in a case in 
which the action is not proper. The set of diagonal matrices form a slice at 
any diagonal matrix, on which all isotropy subgroups are equal, so assumption 
(ii) of the theorem is automatically satisfied. The Ga-^-invariant splitting of Q 
(assumption (i)) is provided by splitting into diagonal and off-diagonal parts. 
Hence from Theorem 8, Eq. (15) spans all foliate vector fields in this case. We 
can think of the foliate vector fields as those for which the evolution of the 
eigenvalues of L depends only on those eigenvalues, and not on the individual 
entries of L. 

More explicit examples of invariant vector fields are given by 

f{L)=p{L)g{trL,trL',...,trL^), 

where p is a real analytic function (which, when extended to matrices, maps 
Q into 0)and g:W M. 

To construct foliate integrators, following Eq. (12) we represent the solution 
as L{t) = U{t)F{t)U{t)-^, U{t) e Gi, F{t) e g, and work with the system on 
Gi X 0, 

if = A{UFU-^)U, f/(0) = /, 

F = f{F), F(0)=L(0), 
Since the symmetry of the second equation is linear, it is preserved by any 
linearly covariant method, such as Runge-Kutta. Runge-Kutta-Munthe-Kaas 
applied to the whole system provides a foliate integrator. The simplest example 
is the following first-order "Lie-Euler" method. Noting that we can take C/„ = 
I, Fn = Ln at each time step (only calculating the updates which move L 
around on the leaves), we get the method 

Un+i = exp(rA(F,)), = F„ + r/(F,), 
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or, in terms of L, 

Ln+i = exp{TA{Ln)){Ln + T f {Ln)) exp(-rA(L„)). 

Example 11 Let M = MT'^'' and let the matrix Lie group G C GL{n) act on 
M by left multiplication, A([/, A) = UAioiU E G, A E M. As in the previous 
example, in the case that the action is proper (e.g. when G is compact) we 
can write the foliate vector fields in the form 

A^g{A)A + f{A) 

for some functions g: M q and f: A ^ A, where / is invariant, i.e. f{UA) — 
Uf{A) for all U E G. As before it is difficult to find all invariant vector fields. 
Examples are f{A) = AV{A^A) for any V:Rp''p Rp^'p, for G = SO{n), 
and f{A) = AV{dctA) for any V:R ^ Rp^'p for G = SL{n). The planar 
systems of Eqs. (1,2) are examples with n = 2, p = 1, and G = 5*0(2). For 
G = SO{n), the leaves are the sets A'^A = const., one of which is the Stiefel 
manifold A^A = I. 

To construct foliate integrators, following Eq. (12) we represent the solution 
as A{t) = U{t)F{t), U(t) e G, F{t) e BI'^'p, and work with the system 

U = g{UF)U, U{0) = I, 
F = f{F), F(0) = A(0), 

As in the previous example, RKMK applied to the whole system provides a 
foliate integrator. The Lie-Euler method is 

An+l = eMM^n))(An + t/(A„)) (16) 

An example of this method is shown in Figure 2. 



5 Conclusions 

The situation considered here, of vector fields preserving a given foliation, can 
be extended in several ways. First, a system may preserve several different foli- 
ations. This is formalized in the 'multifoliate' structure introduced by Kodaira 
and Spencer [12]: the distributions form a lattice, closed under intersections 
and (an appropriate) join. A simple example on is provided by the two 
distributions S- and for which the multifoliate vector fields have the form 

ax ay ' 
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-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

Fig. 2. Foliate vs. nonfoliate integrators. The ODE of Eq. (1) is integrated by the 
(fohatc) Lic-Euler method Eq. (16) (top), and the (nonfohate) Euler method (bot- 
tom). The 20 initial conditions lie on a circle of radius 2, and four time steps of 0.1 
are shown. In the nonfoliate integrator, the final values do not lie on the reference 
circle shown. Note that the two methods coincide on x = 0, where the tangential 
component vanishes. 

A more sophisticated case occurs in Hamiltonian systems with symmetry, 
which preserve the fohations into group orbits and into momentum level sets, 
and in Poisson and conformal Poisson systems, which preserve in addition the 
foliation into symplectic leaves [18]. We plan to study such systems more in 
the future [19]. 

Second, a foliate system may have extra structure, corresponding to the dif- 
ferent infinite dimensional Lie subalgebras of the Lie algebra of foliate vector 
fields. These so-called nonprimitive Lie subalgebras have not been classified. 
However, many examples can be constructed by considering (i) the vector field 
to lie in some other Lie algebra, as of Hamiltonian or volume-preserving vector 
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fields; (ii) the reduced system on the space of leaves to lie in some other Lie 
algebra (in this case the fohated phase space has transverse structure); and 
(iii) the reconstruction system, considered as a nonautonomous vector field on 
a leaf, to lie in some other Lie algebra. (It is the relationships between the 
reduced and reconstruction systems which have not been classified.) Cartan 
[5] provides an interesting list of all infinite-dimensional Lie pseudogroups in 
two dimensions. For example, in the group 



the foliation y — const, is preserved, and the reconstruction dynamics lie in 
the fractional linear group. But in the group 



what structure is preserved? In all of these cases it makes sense to ask how the 
structure affects the dynamics (e.g., what are its homeomorphism invariants), 
how it can be detected in a given system, and how integrators that lie in the 
corresponding group of diffeomorphisms can be constructed. Prom this point 
of view the results presented here are just a beginning. 
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